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THE DISCRETE ANALOGUE OF THE OPERATOR — AND 

ITS PROPERTIES 
Kh.M.Shadimetov 

Abstract 

In this paper the discrete analogue D m [j3] of the differential operator d 2m /dx 2m is constructed 
O^l and its some new properties are proved. 
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1 Main results. 



First S.L.Sobolev [1] studied construction and investigated properties of the operator D^[(3], which 

is inverse of the convolution operator with function = h n G m (hH j3). The function D^{(3] of 

discrete variable, satisfying the equality 

h n D$\0\*G$\0\ = 5\0\ 



is called by the discrete analogue of the polyharmonic operator A m . S.L.Sobolev suggested an algo- 
rithm for finding function D\J[p\ and proved several properties of this function. In one dimensional 

IT) . d 2m 

case, i.e. the discrete analogue of the operator —— = — was constructed by Z.Zh.Zhamalov [2, 31. But 

dx Zm 

there the form of this function was written with m + 1 unknown coefficients. In works [4, 5] these 

O d 2m 

coefficients were found, hereunder the discrete analogue of the operator was constructed com- 

dx zm 

pletely. 

In this paper we give the results of works [4-7], concerning to construction of the discrete analogue 

k>( , ^2m 



D m [f3] of the operator — — = — , and discovery of its properties, which early were not known. 

dx Zm 

Following statements are valid. 

d 2m 

Theorem 1. The discrete analogue of the differential operator — 2m have following form 

m_1 (^ \ ^2m+l 
1 + E e \ X S for |/3| = 1, (1) 

where E a (X) is the Euler polynomial of degree a, Xk are the roots of the Euler polynomial E2 m -2{X), 
in module less than unity, i.e. \Xk\ <1, h is the step of the lattice. 



m[Pi ~ h 2m 
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Property 1. The discrete analogue D m [(3] of the differential operator of order 2m have represen- 
tation 

, -iM m ~ l \l/3|+m-2 

n k=1 ^2m-2\ A k) 

where A^f/S] = (_i) fe + m ^ ^ k/^ ~ ^ ^ s ^ mme ^ c difference of order 2m. 



k=—m 

Property 2. The operator D m {(3] and monomials [(3} k = (h(3) k are connected as 

for 0<*;<2m-l, 



D m \P\ \P\ k - | (2m) , for k = 



'0 for 2m + 1 < fe < 4m - 1, 



(2m)! 

Property 3. T/ie operator D m [j3\ and the function ex.p(2irihp/3) connected as 

J2D m [/3] exp(2mhpp) = 



\mr)2mfr) i \ I t— 2m „• 2m 



(-l) m 2 im (2m-l)!r 2m sin M (i/tp) 

m-~2 

E 

fc=0 



2 £ al 2m - 2) cos27rMm - 1 - *0 + a^ 2) 
w/iere 

a^=B-l^( 2 7)(fc + l- J -)*»- 1 

3=0 V J 7 

are £/te coefficients of the Euler polynomial E2 m -2(X). 

2 Lemmas. 

As known, Euler polynomials E k (X) have following form 



\E k (X) = (l-X) k+2 D k ^^, (4) 



(2) 



^(4m)!* 2 m for , = 4m _ (3) 



where 

D = \-jr, D k = X^D k -\ 
d\ dX 

In [8] was shown, that all roots of the Euler polynomial E k {X) are real, negative and different: 

X[ k) < X (k) < ... < X {k) < 0. (5) 
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Furthermore, the roots, equal standing from the ends of the chain (5) mutually inverse: 

A? } • A<2 W = 1. (6) 

k 

If we denote E k {X) = ^ a^X s , then the coefficients of Euler polynomials, as this was shown by 

s=0 

Euler himself, are expressed by formula 



4 fc) = E(-i) J (^ 2 )( s +i-i) fc+1 - 

3=0 V ' 



From the definition E k (X) follow following statement. 

Lemma 1. For polynomial E k (X) following recurrence relation is valid 



E k (\) = (kX + l)£ fc _i(A) + A(l - X)E' k _ 1 (X), (7) 



where Eq(X) = 1, fe = 1, 2, .... 

Lemma 2. XTie polynomial E k (X) satisfies the identity 



E k (X) = X k E k , 



or otherwise = a k k \, s = 0,1, 2, k. 

Proof of lemma 1. From (4) we can see, that 



(8) 



tf fc _i(A) = X-Hl ~ X ) k+lL,k ~ 1 ^4xj2- W 
Differentiating by A the polynomial E k -i(X), we get 

4-i(A) = "(1 " A) fc A" 2 (fcA + l)^" 1 + AAy. 

Hence and from (9) we obtain, that 

(kX + l)£ fc -i(A) + A(l - X)E' k _ 1 (X) = (kX + lJA-^l - A ) fc+1 ^ _1 7TrX)2- 

_(1 _ A) fe+1 A _1 (A;A + l^^A-^ + E fc (A) = £7 fc (A). 

So, lemma 1 is proved. 

Proof of lemma 2. The lemma we will prove by induction method. When k = 1 from (4) we 
find 

E 1 (X) = A + 1. 

We suppose, that when k > 1 the equality al fc ^ = aj^,^ , n = 0, 1, A; — 1 is fulfilled. We assume, 
(k—l) 

that a n = for n < and n > k — 1 . 
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From (7) we have 

4 fe ) = (s + l)^*" 1 ) + (k - s + l)^, 
then, using assumption of induction, we get 

41 = (*-- + l)4- s 1} + (- + lja^i = (* - * + ljafc 1 ' + (a + lja?" 1 ) = a<*\ 
and lemma 2 is proved. 

3 Proof of theorem 1 

For this we will use function 

9m — 1 

x signx 
° m{x} = 2-(2m-l)!' 

To this function we correspond following function of discrete argument: 

Um W- 2-(2m-l)! • 
Here we must find such function D m [/?] , which satisfies the equality 

hD m \fi]*G m \p] = S\p]. (10) 

According to the theory of periodic generalized functions and Fourier transformation in them 
instead of function D m [(3] it is convenient to search harrow shaped function [1] 



The equality (10) in the class of harrow shaped functions goes to equation 

h U m (X)* Urn (x) = 5(x), (11) 

where 

G m { x 



It is known [1] , that the class of harrow shaped functions and the class of functions of discrete vari- 
ables are isomorphic. So instead of function of discrete argument D m [f3] it is sufficiently to investigate 

the function D m {x), defining from equation (11). 

Later on we need following well known formulas of Fourier transformation: 



F [f(p)] = J f(x)exp(2iripx)dx, 
F^ifip)] = J f(x)exp(-2iripx)dx, 
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F[f(x)*<p(x)]=F[f(x)]-F[p(x)], 
F[5(x)] = 1. 

Applying to both parts of (11) Fourier transformation, we get 

F[D m (x)} ■ F[h ¥ m (x)] = 1. 
Fourier transform of h G m (p) is well known periodic function, given in R with period h^ 1 

This formula is obtained from the equalities 



(12) 



(13) 



and 



F[GM] = j2^W^ ([1 ' P - 729]) 
¥ m (x) = G m (x)^(5(x-/i/3). 



Hence, taking into account (12), we get 



F[D m (p)] 



(2vr) 2m 



y L 



2m 



-1 



(14) 



The main properties of this function in multidimensional case, appearing in construction of discrete 
analogue of the polyharmonic operator, were investigated in [1]. 
We give some of them, which we will use later on. 

1. Zeros of the function F[D m (p)] are the points p = h~ l ft. 

2. The function F[D m (p)} is periodic with period h~ x , real and analytic for all real p. 
The function F[ D m (p)} can be represented in the form of Fourier series 



F[D m (p)] = ^2D m [(3]exp(27rihf3p), 



(15) 



where 



D m [f3]= / F[D m (p)]exp(-2irih(3p)dp. 



Applying inverse Fourier transformation to the equality (15), we get harrow shaped function 

D m (x) 



(16) 



(17) 
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Thus, D m [j3] is searching function D m [(3] of discrete argument or discrete analogue of the operator 
d 2m 

For finding the function D m [j3] calculation of the integral (16) inadvisable. We will find it by 



dx 2m 
following way. 

By virtue of known formula 



y— 

t-^ (rt — i 



IT 2 



(p-/3) 2 sin 2 7rp 
and from the formula (13) we get 

T—T —1 1 

F[h Gi(p)] = 7^E fa-Hffla = 



(2tt) 2 V (P - W 4 sin2 *P h ' 



Hence by differentiating we have 

d_ 
dp 



Thus continuing further, we obtain 

d 2m ~ 2 — (2m - 1)! ^ 1 

dpP^ 2 [ G 1 ~ (2^) 2m ^ (p - /i^ 1 /?) 2 ™ 

= (-l) m ^(2m - l)!(2^) 2m " 2 F[/i cT m (p)]. 

So, 

(-l) m /i 2 d 2m ~ 2 



^[/i G m (p)] 



2 2m vr 2m - 2 (2m - 1)! dp 2m " 2 Vsin 2 7r/ip 
Consider, the expression 

d 2m-2 



dp 2m 2 \sm 2 irhp 
Using 

exp(irihp) — exp(-irihp) 

sin 7mp = , 

2t 

we have 

d 2m ~ 2 / -4 \ d 2m ~ 2 ( exp(2vrmp) 



_4_ 



dp 2m ~ 2 \ v (exp(7Ti/ip) — exp(— irihp)) 2 J dp 2m ~ 2 \ (exp(27ri/ip) — l) 2 / ' 

We will do change of variables A = exp(27ri/ip), then in view of that 

d dX d d d 

~a~ = and T = 27r«/iA— , 

dp dp aA ap aA 

we get 

j2m-2 

- (27rih) 2m ~ 2 D 2m ' 2 , 



dp 



,2m-2 
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where 

D = A—, D 2m ~ 2 = X—D 2m ~ 3 . 
dX dX 



Thus, 



F[h ¥ m (p)] = t^t,,^ 2 ™- 2 A 



(2m -1)! (1-A) 2 ' 
Hence in virtue of (4) we have 



F[/l Gm(p)]= (2m-1)! (1-A)^ - (18) 



From (18), according to (12), we obtain 



F\D (p)}- (2m - 1)] {l ~ X?m (19) 

*[D m [p)\- h2m X E 2m _ 2 (xy U9j 



Now in order to obtain Fourier-series expansion, we will do following. 
We divide the polynomial (1 — A) 2m to the polynomial A£ , 2m-2(A): 



— A - 2m - afc^ + ^, (20) 



2 ^ 2 (2m-2), s 2m " 3 A£ 2m _ 2 (A)' 

A 2^ fls A s 

s=0 

where -P2m-2(A) is a polynomial of degree 2m — 2. It is not difficult to see, that the rational fraction 
2m 2 ( ) . g p r0 p er f rac ti n, i.e. degree of the polynomial P2m-2(A) is less than degree of the 

A£'2m-2(A) 

polynomial XE 2m - 2 (X). Since the roots of the polynomial E 2m _ 2 (X) are real and different, then the 

rational fraction 2m 2 ^ is expanded to the sum of elementary fractions. Searching expansion has 
XE 2m - 2 (X) 

following form 

r> / \ \ A tci—\ . m— 1 . 

^ 2m - 2 {X) = Ap A l)k \ - A 2)k 

XE 2m . 2 (X) A ^ A - Ai, fc ^ A - A 2 , fc ' 1 ' 

where Aq, A\^, A 2 ^ are unknown coefficients, Ai^ are the roots of the polynomial E 2m _ 2 (X), in 
modulus less than unity, and X 2)k are the roots of the polynomial E 2m - 2 (X), in modulus greater than 
unity. By (21) the equality (20) takes the form 



X £ af"- 2) A* 

s=0 



m-1 . m-1 . 

fc=i i,K fc=i Z ' K 



Reducing to the common denominator and omitting it, we get 

(1 _ A) 2m = X 2 E 2m „ 2 {X) - A(2m + a i ^ ) )E 2m ^ 2 (X)+ 



7 



_l4 jp i A i,k^E 2m -2W \^ A 2jk XE 2m -2(X) 

+A E 2m _ 2 (X) + ^ A _ Alfc + 1, A-A 2fc • (23) 
fe=i i,/c fe=i A/c 

Assuming in the equality (23) consequently A = 0, A = Ai^ and A = \ 2j k, we find 

1 = E 2m - 2 (0)A ; (1 - Ai 5 fc) 2m = Ai i fc£ , 2 m ,_2(^i,fc) j 4i,fc; 

Hence 

4 _i 4 (1-Ai, fc ) 2m , (l~A 2 , fc ) 2m 

> ^1^— X 771/ /X V ^2,fc — . 771/ /X \- 

Al,fe^2m-2(. A l,fcJ M,kti 2m - 2 {A2,k) 

Using (6), we have 

fl-T^-) 2 ™ r\ l^m 

/ 1 \ \2m— 1 771/ i 1 \ ' 
A 2,fe- C/ 2m-2vA 2 , fe ' 1 A 2,fe -^2/71-2 V A 2jfe J 

In virtue of (7) we obtain 

T ( 1 - ~x ) E 2m _ 2 (- ) = E 2m -i(- ), 

M,k V A 2,fc/ A 2,fc A 2,fc 

A2,fc(l - A 2 ,fc)-E , 2m-2(A2,fc) = ^2m-l(A2,fc), 

hence 

^2771-2 U J — 



^2m-2(^2,fc) 



A2,fc A2,fc — 1 

^2m-l(A2,fe) 



A2,fc(l - A 2 ,fc) 
From here application of the lemma 2 gives 



4 4 _ (1 ~ Al, fc ) 2 " t+1 f24 x 



Since | Ax,fc | < 1 and |A 2; fc| > 1, then 



m-l , m-l - 



^A-Ai fc Z -^A-A 2 fc 

k=l ' fc=l ' 

can be represented as Laurent series on the circle | A 2 1 = 1 : 

m— 1 < .. m— 1 « 1 m— 1 oo ✓ . /3 



k=l L ' K k=l 1 A k=l 0=0 v 7 



m-l . m-l , m-l . 00 / \ \ Z 3 

^2,fc = A 2) k = \ - ^2,fc / A \ 
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Putting (25), (26) to (22) and taking into account A = exp(27ri/ip) from (19), (20), we obtain 

(2m - 1)! 



F[D m (p)\ = 



m—l 



exp(27ri/ij>) — 2m — <4^_ 3 2 ^+ 



+ exp(-2irihp) + A lt k^2 Xl k exp(-2nihp(f3 + 1))- 

k=l V 0=0 

-^2,fcA-j^ exp(2irih P p) 



0=0 



Thus, searching Fourier series for F[ D m (p)] have following form 

F[D m (p)] =^2D m \p]exp(2mhpp), 



where 



D m [f3] 



(2m- 1)! 



m—l 

k=l 

m—l 

i + £ M, k 

k=l 

m—l 

22m— 1 

k=l 
m—l 

1 ~ £ A 2,kKk 
k=l 
m—l 

~ £ ^2,feA^ 
fc=l 



1 for /?<-2, 

for /3 = -1, 

£ A 2 ,k^i for = o. 



0-i 



for (9 = 1, 
for p > 2. 



With the help (24) the function D m [/3] we rewrite in the form 



(2m - 1)! 

P™ 



m—l /-i \ \2m+l \ 
^ Ajfc-E2m-l(Al,fe) 

(1 - Ai )fc ) 



We note, that 

Theorem 1 is proved completely. 



fc=i 

_22m— 
A, 



-E , 2m-l(Al j fe) 

m ~! /i \ \2m+l 
1 + (1 ~ M,k) 

~_ Al,jfc-E , 2m-l(Al,fe) 



for > 2, 
for |/3| = 1, 
for = 0. 
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4 Proofs of properties. 

Proof of property 1. 

Following is takes placed 



^[A^ (p)] = -4sin 2 7rp/i. 



Indeed, 



A 2 (x) = 6 (x + l)-2 5 (x) + 2 S {x - 1) = S(x + h) - 28{x) + 8{x - h). 
By definition of Fourier transformation we have 

F[/S.2 (p)] = J ex.p(2iTipx) A 2 {x)dx = — 4 sin 2 irph. 
Hence consequently we obtain 

m times 

/ " V 

T—rlm] r—rll] ,— ,[1] ,— ,[1] 

F[A 2 (p)]=F[A 2 (p)* A 2 (p)*...*A 2 (p)] = (-4rsin 2m 7r^. 



Immediately we have 



(1 ~ A) ; 
-4A 



sin 2m irhp. 



By virtue of (27) and (28) the formula (19) takes form 

, v, (2m -1)! A™" 1 

Now expanding rational fraction -g^ — to the sum of elementary fractions, we have 



(27) 
(28) 

(29) 



A 



m— 1 



-E'2m-2(A) 



TO— 1 

E 

fc=i 



B 



l,k 



B 



2,k 



A- A 



where 



B 



l,k 



A l,fc 



E 2m-l{^l-,k)' 



B 



2,k 



A-A 2 



\TO — 1 
\k 

E 2m-2( X 2,k) 



(30) 



Since | Ax,jfc I < 1 an d |A 2 ,fc| > 1> then expanding x Bl ^ - and A f - to the Laurent series on the circle 
|A| = 1, we find 



B 



l,k 



l.k 



A - Ai )fc A i _ dhh. 



E 



0=0 



Ai,fc 
A 



B 



2, A; 



B 



2,fe 



B- 



A-A 2 . fc A 2)fc (l-^; 



2,k 
^2,k 



E 

0=0 



A 
A 2 j 



(31) 
(32) 
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On the strength of (6), (7), (8), (31), (32) the equality (30) takes the form 

X m-i 

^2m-2(A) ^ E 2m-2( X k, ^ 

So, from (29) we get 



m-l \m-2 



no m w] = M E E Af'A^i a r wi. 

n k=l ^2m-2^k) p 

Applying to (33) inverse Fourier transformation, after some simplifications we obtain 



(2m -1)!^. 



h 2 



m-l ,|/3|+m-2 



(33) 



(34) 



According to definition of harrow shaped functions from (34) we get the statement of property 1. 

Proof of property 2. 

The equality (2) proved in [1]. 
Here we will prove (3). 

From (15), (17) we obtain 



F[D m (p)] =J2 D m[P]exp(2Trih/3p). 



(35) 



Using expansion of exp(27ri/ip/3) and by formula (2) in cases < k < 2m — 1 and k = 2m, we will have 

{2mhpp) k 



F[D m ( P )]=Y, D mW]Y,- 



k=0 



k\ 



/3 k=2m+l 

On the other hand, from (14) we get 

F[D m (p)] = {2vi) 2m 
Thus, on the strength of (36) and (37) we have 



E 



1 R\2m 



(36) 



(37) 



E D ™1® £ (2m tf f + (2^) 2m = (2^ 



fc=2m+l 



fc! 



E 



{p-h- l (3) 2m 
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or 



£a»[/3] £ 



(2irihp(3) k 



fc=2m+l 



fc! 



Consider 



ip(h,p, m) = 1 



= -(2vrip) 2m | 1 - 
(p/i) 2m 



E 



2m 



(p/i - 7) 2m 



E 



1 - 



i + ( P /o 2m E 



(p/i - 7 ) 2m 



+ 



-i 



^[(p/i- 7 ) 2m (p/i + 7) 



1 2m 



-1 



= 1 - 



7=1 

Hence choosing p such, that 



■2m 



(1 _ PJl^-Zm _)_ (1 _)_ P -j-2rn 



7 



7 



7=1 



-2m 



(X — PJl^- 2m _|_ (X + ^)~ 2m 



7 



7 



< 1, 



expanding the fraction — to the series of geometric progression, we have 



1 + Q 

oo 

ij{h,p,m) = (ph) 2m Y J 



-2m 



7=1 



(X _ ^j- 2 *" + (1 + ^)~ 2m 



7 



7 



+ 



+0((M 4m ) = 2(p/i) 2m ]T 7 _2m + 0(/i 2m+1 )- 

7=1 

The left part of (38) we rewrite in the form 



\k 00 



P k=2m+l " n=l 



sr^ (2TripB) 2m+n 
where fj,(p, n, m) = 22 h2m °^\ ^m + ny. d ° eS n0t depend on /l ' 

since h 2m D m [j3] does not depend 

h, that clear from (1). 

Comparing right hand sides of (39) and (40), from (38) we have /i(p, n, m) = for n = 1, 2, 2m- 
1, i.e. for k = 2m + 1, 4m — 1 



(38) 



(39) 



(40) 



n(p,2m,m) = Y,h 2m D m [P] 



(27rip(3) 
(4m)! 



4m 



(-l) m+1 2(2 ff ) 2 V m E ^ _2m - 

7=1 
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which proves the property 2. 



Proof of property 3. From (19) and (35) we have 




2m 



^Dmtf] exp{2mhpp) = 



Using (8), (40) and A = exp(27rz/ip/3), after some simplifications we get property 3. 
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